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Abstract 

The generator-coordinate method is a flexible and powerful reformulation of the variational 
principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation 
of density-functional theory, excitation energies can be obtained from ground-state density func- 
tional. As a viability test, the method is applied to ground-state energies and various types 
of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results 
are compared to a variety of alternative DFT-based approaches to excited states, in particular 
time-dependent density- functional theory with exact and approximate potentials. 
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I. INTRODUCTION 

Density-functional theory (DFT) l[ [2], [J is the most widely used first-principles method 
for the calculation of ground-state properties of a wide variety of atomic, molecular and solid 
systems. Common local and nonlocal density functionals make highly precise predictions 
for ground-state densities, energies and quantities that can be derived from these. Efficient 
algorithms for solving the Kohn-Sham (KS) equations are implemented in many common 
electronic-structure codes, and provide computational access to ground-state densities by 
means of an auxiliary set of single-particle orbitals. 

The Kohn-Sham equations, as well as most common density functionals, were designed 
with ground-state properties in mind. The underlying theorems of density-functional theory, 
however, are more general. The usual proofs of the Hohenberg-Kohn theorem by contra- 
diction fl or by constrained search Q guarantee that the ground-state wave function is a 
functional of the ground-state density, but it takes only one simple additional step to prove 
that the density also determines the external potential and thus all wave functions 

and energies, including those of excited states s|. 

DFT thus holds the promise to become a versatile and powerful tool for the calculation 
of excited-state energies, too. In practice, however, this promise turned out to be much 
harder to fulfill than for the ground state. In fact, the most commonly employed DFT 
approaches to excited states are formulated in the conceptually rather different frameworks 
of time-dependent DFT [9| or ensemble DFT [10]. For a comparison of these methods, and 
a variety of other DFT-based approaches to excited states, see Ref . 11[ . What TDDFT and 
ensemble DFT have in common is that by design they go, from the outset, beyond static 
(ground-state) DFT. In this paper we propose and test an alternative method that allows 
one to extract excited-state energies from ordinary ground-state functionals. 



II. THE GENERATOR COORDINATE VARIATIONAL PRINCIPLE 

The key ingredient in this development is the generator coordinate variational principle. 
The generator-coordinate method (GCM) arose in nuclear physics, as a way to build collec- 
tive behaviour of nuclei into a trial wave function written in terms of single-particle orbitals 
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12j. In this method, the (nuclear) many-body wave function \1/ is cast as 

fy(x ll ..,x N )= daf(a)$(a;xi,..,x N ), (2.1) 



where the $ are auxiliary wave functions arising from a deformed Hamiltonian, and the 
degree of deformation is characterized by the deformation parameter a, which plays the role 
of a generator coordinate. It is a crucial feature of the method that the generator coordinate 
appears in the auxiliary (seed) functions $, but not in the wave function Variation of 
the energy with respect to the weight function /(a) leads to an integral equation whose 
eigenvalues are the many-body energies of the system, 

J da' [K(a, a') - ES(a, a')) f(a') = 0. (2.2) 

Equation (12.21) . known as the Griffin-Hill- Wheeler (GHW) equation, contains Hamiltonian 
and overlap kernels K (a, a') = (&(a)\H\$(a')) and S(a, a') = ($(a)| $(«')) of standard 
form, and can be solved straightforwardly by discretizing the a integral and obtaining the 
eigenvalues E and eigenfunctions f{a) by matrix algebra. 

The GHW method is also used in quantum chemistry to construct highly precise basis sets 
for Hartree-Fock and Dirac-Fock calculations [13J. In these applications the single-particle 
orbitals are written in GCM form as 

(p(x) = J da f(a)£(a; x), (2.3) 

where is a set of suitable single-body functions, a is identified with a basis set exponent, 
and <p(x) is the Hartree-Fock (or Dirac-Fock) orbital. The resulting GHW equation produces 
numerically defined basis sets of high accuracy 

In 2003 one of the present authors proposed a way to combine the GCM variational 
principle with DFT [ijj], by identifying the deformation potential of the original GCM with 
the KS potential of DFT, the generating function with the KS Slater determinant, and the 
generator coordinate with any parameter in the effective single-particle potential. Hence, 

V(x u ..,x N ) = J da f(a)$(a;x 1 ,..,x N ), (2.4) 

where ^ now is the electronic many-body wave function, and $ is the KS Slater determinant. 



Conceptually, this is much more similar to the original use of GCM in nuclear physics 12] 



than to its more recent applications in quantum chemistry 
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A first implementation of GCM-DFT 14j showed that even simple approximations to 
the full formalism result in ground-state energies that are comparable to those obtained 
from more sophisticated density functionals. Moreover, unlike standard DFT methods, 
GCM-DFT also provides an approximation to the many-body wave function, not just to the 
energies. The final expressions obtained for these wave functions are formally similar to CI 
expansions, but whereas in CI each determinant represents a different excited-state config- 
uration of one fixed Hamiltonian, in the GCM-DFT method of Ref. [ijj] each determinant 
comes from a ground-state calculation of a differently deformed Hamiltonian, and can thus 
be viewed as a resummation of many CI determinants. As these GCM determinants come 
from different Hamiltonians they are not necessarily orthogonal, in contrast with CI deter- 
minants. While at first sight this increases the computational effort involved in constructing 
the GCM-DFT expansion relative to that for constructing the CI expansion, the fact that 
each GCM determinant can be viewed as a resummation of many CI determinants suggests 
that much smaller expansions may be sufficient with GCM than with CI, in particularly if 
the seed functions are cleverly chosen. 

Ref. jl4] also showed that approximations to the energies of excited states could be 
obtained as by-products of the ground-state calculation, but the numerical example given 
there showed that this method performs rather poorly. Here we present a reformulation of 
GCM-DFT which performs much better for excited states - in fact, even in its simplest form 
it is already competitive with standard TDDFT based methods. 



III. GCM-DFT FOR EXCITED STATES 



To introduce the key idea, let us first briefly review the way ground-state energies are ob- 
tained from GCM-DFT [14!]. The seed wave functions (or generator functions) <£>(a; xi, ..,xn) 
are chosen to be the Slater determinants arising from Kohn-Sham orbitals of the N lowest- 
lying levels of the deformed KS Hamiltonian, i.e., are KS A-particle ground states. The 
particular deformation chosen in Ref. [ijj] was 

v a (r) = v ext (v) + v h (t) + av^ DA (r), (3.5) 

i.e., the exchange-only LDA was modulated by a, for which the five-point mesh 
{0,0.5,1,1.5,2} was used. We stress that this is not a Xa calculation, although the last 
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term in Eq. (13.51) is of Xa form: a here is a generator coordinate, which appears as integra- 
tion variable in the integral-representation of the wave function (12.41) . Upon discretization 
of the integral, a assumes more than one value, none of which is fitted to experiment. For 
each a we selfconsistently solve the KS equations, construct $(a) from the resulting orbitals 
and calculate the kernels K(a,a') and S(a,a'). 

The benchmark He ground-state energy E^ e = —2. 904a. u. was reproduced in this way to 
within 1.1%. A further improvement on the ground-state energy, reducing the deviation to 
—0.24%, is reported below. (For comparison, the LDA, when used in the standard DFT way 
without generator coordinates, predicts jr;^ e ' LDA = — 2. 8348a. -u., which deviates by -2.4% 
from Eq £ .) 

In principle, the five eigenvalues arising from the 5x5 matrix problem allow one to obtain 
four excited-state energies, in addition to the ground state, but the lowest-lying excitation 
energy was found [3| to deviate by 16.7% from the exact value, which is too much to be 
useful in practice. 

However, we note that the trial function \l/ inherits its symmetries and structure from 
the seed functions $(a). We can thus target a particular excited state of the many-body 
system by using seed functions of the corresponding excited KS state, and solving the GHW 
equation with kernels obtained from these excited-state determinants or configuration-state 
functions. This procedure can be used for all many-body excitations that have a counterpart 
in the noninteracting KS system. (Excitations not having such a counterpart can be dealt 
with by using a generator coordinate in TDDFT [15j, or by using more sophisticated seed 
functions in static GCM-DFT.) 

Determinants corresponding to excited states of the noninteracting KS system are rarely 
used in DFT, because the KS formalism is tailored to provide ground-state properties. 
However, in the context of GCM-DFT, the KS equation is not used to obtain the ground- 
state density, but to obtain a set of continuously parametrized iV-particle determinants that 
are employed as seed functions <£>(«) in the GCM ansatz Eq. (12. lft . and for this purpose the 
use of excited-state KS determinants is perfectly legitimate. 
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IV. APPLICATIONS TO REPRESENTATIVE TWO- AND THREE-ELECTRON 
SYSTEMS 



In this section we report numerical results for different types of excitation energies (sin- 
glet/doublet/triplet) as well as ground-state energies, for various two- and three-electron 
systems from the He and Li isoelectronic series, compared to a variety of other computa- 
tional approaches, including TDDFT. 



A. Excited-state energies for atoms from the He isoelectronic series 

First, we present illustrative applications to excited states of the Helium isoelectronic 
series, and to the lowest singlet and triplet excitations of the He atom. We chose the 
deformation potential to be 

v a (v) = v ext {v) + v H (r) + a[v L x DA {v) + v L c DA {v)l (4.6) 

i.e., let the generator coordinate modulate the LDA for exchange and correlation. Many 
other choices are possible, e.g., introducing a in the correlation potential, Hartree potential, 
external potential, kinetic energy term, or in the angular or spin-dependent part of the 
orbitals, or using other functionals than LDA. Each of the resulting deformations has its 
own physical meaning and consequences. Each also implies a, generally distinct, range of 
values of a, and, consequently, different meshes for discretizing the integral equation. 

Some choices of where to place a in the KS equations, and some discretization schemes 
for the a integral, give excellent results for specific quantities or some particular system. 
Below we do not report such special choices or best results, but focuse on a simple and 



generally usable scheme [161 ] . Systematic exploration of the many possibilities that arise 
upon combining the GCM idea with DFT remains work for the future. 

In Table (TJ we report lowest triplet excited-state energies of the Helium isoelectronic 
series, obtained from solving the discretized GHW eigenvalue equation with v a (r) chosen 
as in Eq. (I4.6p . on the mesh {4.5,5,5.5,6,6.5}. (Denser meshes lead to only marginal im- 
provements of the results, or may even worsen them if the resulting overlap matrix S(a, a') 
becomes singular.) The seed wave functions were constructed from the KS orbitals arising 
selfconsistently in potential (14. 6p . selected to form the lowest KS excited state of triplet sym- 
metry. Near-exact theoretical data are also reported, and show that the resulting excited- 



state energies are surprisingly close. Table [TT] makes the same comparison for the lowest 
singlet excitation, 2 1 5'. 

For spectroscopy, the interesting quantities are not primarily the excited-state energies, 
but the excitation energies. In Table II I II we thus present selected excitation energies of the 
He atom: the lowest singlet and triplet excitation energies, and the singlet-triplet splitting. 
GCM-DFT can be used as a stand-alone method to obtain excitation energies if it is used to 
calculate the energies of the ground state and of the excited state of interest. Alternatively, 
it can be used as an add-on method by adding excitation energies obtained from GCM-DFT 
to ground-state energies obtained by traditional methods. Both procedures are compared 
in Table III 1 1 to six other computational schemes and near-exact benchmark data. 

The columns KS exact and KS LDA report the KS single-particle gap between the highest 
occupied and lowest unoccupied KS eigenvalue, arising from the exact [l8j KS potential and 
from the LDA potential for He, respectively. Note that neither the LDA nor the exact KS 
eigenvalues predict any singlet-triplet splitting. The column labeled ASCF reports LDA 
total-energy differences between the ground state and the total energy obtained from a 
standard KS calculation fixing the occupation of KS levels at that of the corresponding KS 
excited state. 

The two columns labelled TDDFT**EXX and TDDFT**ALDA report data obtained in 
Ref. is| from TDDFT using adiabatic exact-exchange (AEXX) and adiabatic LDA (ALDA), 
respectively. These TDDFT data were obtained from a numerically exact potential 18], 
approximating only the xc kernel, and in this sense are not fully representative of standard 
TDDFT calculations. ALDA results obtained from the standard TDDFT procedure, using 
approximate potentials and kernels, are reported in the column labelled TDDFT ALDA. 
Our LDA and TDDFT ALDA data were obtained with the GAUSSIAN 03 Q program, 
using the VWN5 parametrization of the^LDA and the aug — ccpV5Z basis set, and with the 
mesh-based atomic DFT code opmks j^. 

The column labeled GCM-DFT* reports differences between GCM-DFT excited-state 
energies and the exact ground-state energy, whereas the column labelled GCM-DFT reports 
differences between GCM-DFT excited-state energies and GCM-DFT ground-state energies, 
the latter being obtained from Eq. (14.61) on the mesh {4.7,5.05,5.4,5.75,6.1}. The He 
ground-state energy obtained on this mesh is — 2. 897a. w., which deviates from the exact 
result -2.904a.u. by -0.24%. The column labeled GCM-DFT* thus measures the add- 
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on performance of generator-coordinate DFT, in our simple implementation, purely for 
excited states, while the column labeled GCM-DFT quantifies its stand-alone performance 
for excitation energies. 

As the data in the lines labelled 3 S and 2 S of Table II II I show, GCM-DFT, both in its 
add-on and its stand-alone version, is capable of producing similar or better singlet and 
triplet excited-state energies than the other tested methods. As the line labelled A shows, 
it also produces realistic singlet-triplet splittings. Still better singlet-triplet splittings are 
obtained from TDDFT/ALDA, both with approximate and exact single-particle potentials, 
but this improvement is due to error cancellation between the energy of the singlet and the 
triplet excited state, each of which individually has larger errors. 

Note that we have not separately optimized the mesh for singlet excitations and triplet 
excitations, but used the same five values of a for both. In other calculations, we have 
obtained better excitation energies by employing different meshes for different states, but 
our purpose here is to keep the procedure as simple and generally applicable as possible, 
and we thus used the same set of a's for all excited states of He. 

B. Ground-state and excitation energies for atoms from the Li isoelectronic series 

As an initial exploration of three-electron systems, we now present GCM-DFT calcula- 
tions of the ground-state energy and excitation energies of atoms from the Li isoelectronic 
series. In principle, one could choose a new placement of the generator coordinate in the KS 
equation, and a new mesh, for each new system. In this initial exploration, we maintained for 
Li the choices made for He. Ground-state energies obtained in this way are reported in the 
first column of Table HVl labelled 'He mesh'. Clearly, the deviations from the experimental 
energies are much larger than they were for the He series. 

Next, in recognition of the fact that the average radius (r) of an atom shrinks as Z 
increases (for one-electron atoms, (r) oc 1/Z), and taking into account that the generator 
coordinate in Eq. (14.61) modifies only the radial wave function, we constructed downscaled 
radial meshes for Li by applying a simple power-law scaling factor to the mesh for He, 



where the exponent 3/4 was chosen after some experimentation (but not optimized varia- 
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(4.7) 



S 



tionally). This scaling approach is quite successful: The second set of data in Table HVl has 
much reduced deviations from the benchmark data. 

Finally, instead of adjusting the exponent of the mesh scaling law (j4.7p . we also con- 
structed a new discretization mesh for neutral Li, which was subsequently used also for the 
ions. The free parameters here are the starting value and the increment, i.e. one more 
than in the scaling approach. Results are shown in the third data set in Table IIV} labelled 
'emp', for empirical (but not fitted). As expected, the agreement with benchmark data is 
further improved. We stress that in using these modified meshes we still employed the same 
generator coordinate as in Eq. H4.6f) and still constrained the wave function to be a sum of 
five terms, corresponding to five values of a, i.e., our choice of mesh is an exploration of 
the capabilities of the GCM method, not an arbitrary construction designed to produce the 
best possible agreement with reference data. 

Turning now from ground states to excited states, Table IVl presents KS, ASCF, TDDFT 
and GCM-DFT data for excitation energies corresponding to the ground-state to doublet- 
excited-state 2 S transition ls^s 1 -> ls 2 3s l . As in TableHVl the first group of Li GCM-DFT 
data represent the stand-alone approach, employing the He mesh, a power-law scaled He 
mesh and a mesh whose initial value and step size were chosen to deliver values close to 
the benchmark data for neutral Li, but subject to the constraints of same placement of the 
generator coordinate and same number of mesh points. The only difference is that for the 
excited states the scaling law exponent 1/3, leading to the mesh 



was found to be more suitable than 3/4. The second group of GCM-DFT data, labelled 
GCM-DFT*, reports corresponding results from the add-on use of GCM, in which only the 
excited-state energy was obtained from GCM-DFT, the ground-state energy employed was 
the exact one. 

As for He, GCM-DFT produces realistic excitation energies, but the agreement with the 
reference values is less good than for He. Nevertheless, since we have not fitted, and neither 
variationally optimized, the mesh of a values, but just explored some simple and generaliz- 
able discretization schemes, the agreement achieved is still rather encouraging. Moreover, 
we note that the smaller errors TDDFT achieved for the 2 2 S — > excitation energies are 
a consequence of error cancellation between the ground-state energy and the excited-state 




(4.8) 
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energy, each of which is less well reproduced by TDDFT/ALDA than by the GCM-DFT 
calculations. This is reminiscent of the error cancellation that occured between singlet and 
triplet excited- state energies for He. 



V. CONCLUSIONS 

Whether GCM-DFT will ever be competitive with more established TDDFT or CI ap- 
proaches to excited states is, at present, an open question, but from the present analysis it 
seems safe to conclude that, as a matter of principle, excitation energies can indeed be ob- 
tained from static ground-state density functionals, without requiring time- or temperature- 
dependent generalizations, as anticipated by the original Hohenberg-Kohn theorem. 

In our opinion, what is missing to turn the presently proposed approach into a viable and 
transferable method for calculating ground-state and excitation energies (and other prop- 
erties 241 ]) of many-electron systems is a systematic and nonempirical way of constructing 
suitable meshes for discretizing the a integral in the GHW equation. Different ways of 
achieving this are conceivable: 

(i) Much denser meshes, employing hundreds or more values of a, would lead to an ever 
more faithful representation of the GHW integral itself, but in this case much care must be 
taken in inverting the overlap matrix, because closely spaced values of a lead to near linear 
dependence of rows and columns and thus to potentially ill-conditioned eigenvalue problems. 

(ii) Variational optimization of a set of a values may be useful for the ground-state energy, 
and produces a mesh that can be used as a starting point also for excited states. 

(iii) Scaling laws, as our present Eqs. (14. 7p and (14. 8p . provide a way to use chemical or 
physical intuition to pre-select a suitable set of values of the generator coordinate without 
data-fitting or optimization. 

This discussion shows that the great strength of the GCM approach is also its Achilles 
heel: the fact that the generator coordinate does not explicitly appear in the wave function 
makes it a particularly powerful and flexible formulation of the variational principle, but 
also implies that there is no direct way to make an a priori reasonable choice of placement 
and range of the generator coordinate. 

We end this paper by noting that, very recently, generator coordinates have also been 
introduced in TDDFT in order to describe retardation effects and memory by means of a 
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time-dependent generalization of the GHW equations [15j]. A further interesting possibility 
is to introduce a generator coordinate also in the TDDFT approach to excited states in order 
to obtain the time-dependent many-body wave function, or improved excitation energies. 

This work was supported by FAPESP, CAPES and CNPq. We thank C. A. Ullrich for 
useful discussions. 
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TABLE I: Energy of the lowest excited triplet state, 2 3 S, of ions from the He isoelectronic series 
obtained by GCM-DFT from the LDA, compared to available near-exact theoretical data [17]. 





rpGCM—DFT 

1 J 


— F, ex 

111 


/U U-C- V IdLlUll 


He 


2.173 


2.175 


-0.092 % 


Li+ 


5.109 


5.104 


0.098 % 


Be 2+ 


9.294 


9.289 


0.054 % 


B 3+ 


14.73 


14.72 


0.068 % 




21.42 


21.41 


0.047% 



TABLE II: Energy of the lowest excited singlet state, 2 1 S, of ions from the He isoelectronic series 



obtained by GCM-DFT from the LDA, compared to available near-exact theoretical data 
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TABLE III: Lowest triplet and singlet excitation energies and singlet-triplet splittings of the He 
atom, obtained by eight different calculational approaches, described in the main text, and near- 
exact benchmark data. Second line for each excitation: percentage deviation from benchmark data 
The last two rows report the singlet-triplet splitting A (multiplied by 100 for legibility) and 
its percentage error. 
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0.7460 0.6218 0.7292 


0.7659 


0.7678 


0.6373 


0.7667 


0.7600 


0.7578 


% 


-1.6 -17.9 


-3.8 


1.1 


1.3 


-16 


-1.2 


0.29 




100 A 





1.46 


4.52 


3.27 


2.73 


3.55 


3.60 


2.93 


% 


-100 -100 


-50 


54 


12 


-6.8 


22 


23 





TABLE IV: Ground-state energy of ions from the Li isoelectronic series obtained by GCM-DFT 
from the LDA and the respective deviation (%), compared to near-exact reference data [2l|]. In 
the first set of data we use the same mesh already employed for the He ground state, in the second 
set we used the scaled He mesh of Eq. (|4.7[) : a = {3.5; 3.9; 4.3; 4.7; 5.1} , as described in the main 
text, and in the third set we employed the empirical mesh a = {3.47; 3.73; 3.99; 4.25; 4.51}. For 
comparison, the LDA prediction for Li is —Eq DA = 7. 3440a. u., which deviates from the exact 
result by -1.79%. 



2 S 


rpGCM—DFT 


rpGCM—DFT 


rpGCM—DFT 


_ E exact £]J 




He mesh 


Eq. gZD 


emp. 




Li 


7.3179 (-2.14%) 


7.4282 (-0.667%) 


7.4742 (-0.0522%) 


7.4781 


Be 1+ 


14.231 (-0.656%) 


14.269 (-0.391%) 


14.220 (-0.733%) 


14.325 




23.146 (-1.19%) 


23.367 (-0.248%) 


23.335 (-0.384%) 


23.425 


C 3 + 


34.699 (-0.221%) 


34.749 (-0.0776%) 


34.681 (-0.273%) 


34.776 
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TABLE V: First two lines: 2 2 S -> 3 2 5 excitation energy of the Li atom, corresponding to the 
ls 2 2s l -► ls 2 ?,s l single-particle transition, obtained from KS, ASCF, TDDFT/ALDA and various 
GCM-DFT approaches, and its percentage error relative to the near-exact value [2]. Second 
set of two lines: Ground-state energy of the Li atom and its percentage error. Third set of two 
lines: Energy of the 3 2 S excited state and its percentage error. In the column labelled KS we 
report, instead of total ground-state and excited-state energies, the KS single-particle energies of 
the 2s and 3s state. The ground-state energy in the columns labelled GCM-DFT*, representing 
the add-on use of GCM-DFT is, by definition of the add-on procedure, the exact one, while in the 



columns labelled GCM-DFT, representing the stand-alone use of GCM-DFT, it was taken from 
the GCM-DFT values of Table ED 



Li 


KS ASCF 


TDDFT 


GCM-DFT 


GCM-DFT 


GCM-DFT 


GCM-DFT* 


GCM-DFT* 


GCM-DFT* 


exact 




LDA LDA 


ALDA 


LDA/He 


LDA/scal. 


LDA/emp. 


LDA/Hc 


LDA/scal 


LDA/emp 




2 2 S -> 3 2 5 


0.1156 0.1199 


0.1140 


0.2670 


0.08080 


0.1240 


0.4272 


0.1307 


0.1279 


0.1240 


% 


-6.77 -3.31 


-8.06 


115 


-34.8 


0.00 


245 


5.40 


3.15 




-E(2 2 S) 


(0.11628) 7.3440 


7.3439 


7.3179 


7.4282 


7.4742 


7.4781 


7.4781 


7.4781 




% 


-1.79 


-1.79 


-2.14 


-0.667 


-0.0522 










-E(3 2 S) 


(0.00066) 7.2241 


7.2299 


7.0509 


7.3474 


7.3502 


7.0509 


7.3474 


7.3502 


7.3539 


% 


-1.77 


-1.69 


-4.12 


-0.0884 


-0.0503 


-4.12 


-0.0884 


-0.0503 
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